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Abstract 

A Lagrangean and a set of Feynman rules are presented for non-relativistic QFT's 
with manifest power counting in the heavy particle velocity v. A regime is identified 
in which energies and momenta are of order Mv. It is neither identical to the ultrasoft 
regime corresponding to radiative processes with energies and momenta of order Mv^ , 
nor to the potential regime with on shell heavy particles and Coulomb binding. In 
this soft regime, massless particles are on shell, and heavy particle propagators become 
static. Examples show that it contributes to one- and two- loop corrections of scattering 
and production amplitudes near threshold. Hence, NRQFT agrees with the results of 
threshold expansion. A simple example also demonstrates the power of dimensional 
regularisation in NRQFT. 



Suggested PACS numbers: 12.38.Bx, 12.39.Hg, 12.39.Jh. 

Suggested Keywords: non-relativistic QCD, effective field theory, threshold expansion, 

dimensional regularisation. 



^Email: hgrie@phys.washington.edu 



1 Introduction 



Velocity power counting in Non-Relativistic Quantum Field Theories (NRQFT) es- 
pecially in NRQCD and NRQED, and identification of the relevant energy and momentum 
regimes has proven more difficult than previously believed. In a recent article, Beneke and 
Smirnov pointed out that the velocity rescaling rules proposed by Luke and Manohar for 
Coulomb interactions 0], and by Grinstein and Rothstein for bremsstrahlung processes |Q, 
as united by Luke and Savage 0, and by Labelle's power counting scheme in time ordered 
perturbation theory [0, do not reproduce the correct behaviour of the two gluon exchange 
contribution to Coulomb scattering between non-relativistic particles near threshold. This 
has cast some doubt whether NRQCD, especially in its dimensionally regularised version |Q, 
can be formulated using a self-consistent low energy Lagrangean. The aim of this article is 
to demonstrate that a Lagrangean establishing explicit velocity power counting exists, and 
to show that this Lagrangean reproduces the results in Ref. [§. 

This article is confined to outlining the ideas to resolve the puzzle, postponing more 
formal arguments, calculations and derivations to a future, longer publication which will 
also deal with gauge theories and exemplary calculations. It is organised as follows: In 
Sect. the relevant regimes of NRQFT are identified. A simple example demonstrates 
the usefulness of dimensional regularisation in enabling explicit velocity power counting. 
Sect, ^proposes the rescaling rules necessary for a Lagrangean with manifest velocity power 
counting. The Feynman rules are given. Simple examples in Sect. ^ establish further the 
necessity of the new, soft regime introduced in Sect. |^. Summary and outlook conclude the 
article, together with an appendix on split dimensional regularisation 

2 Idea of Dimensionally Regularised NRQFT 

For the sake of simplicity, let us - following ^ - deal with the Lagrangean 

C = {d^^K)Hd^^K) - M^^l^n + ^{d^A){d>'A) - 2M^4$rA (2.1) 

of a heavy, complex scalar field with mass M coupled to a massless, real scalar A. The 
coupling constant g has been chosen dimensionless. $r will be referred to as "quark" and 
A as "gluon" in a slight but clarifying abuse of language. In NRQFT, excitations with four- 
momenta bigger than M are integrated out, giving rise to four-point interactions between 
quarks. The first terms of the non-relativistically reduced Lagrangean read 

>Cnrqft = $^(i5o + ^-gc,A)<!> + + C2($t$)' + . . . , (2.2) 

where the non-relativistic quark field is $ = -\/2Me'^^*$R and the coefficients q are to be 
determined by matching relativistic and non-relativistic scattering amplitudes. To lowest 
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order, Ci = 1 and C2 = oa'^ap ■ '^^^ non-relativistic propagators are 



where T = po — M = -^ + ... is the kinetic energy of the quark. 

When a Coulombic bound state of two quarks exists, the two typical energy and mo- 
mentum scales in the non-relativistic system are the bound state energy Mw^ and the rel- 
ative momentum of the two quarks Mv (i.e. the inverse size of the bound state) [^. Here, 
u = /37 -C 1 is the relativistic generalisation of the relative particle velocity. Cuts and poles 
in scattering amplitudes close to threshold stem from bound states and on-shell propagation 
of particles in intermediate states. They give rise to infrared divergences, and in general 
dominate contributions to scattering amplitudes. With the two scales at hand, and energies 
and momenta being of either scale, three regimes are identified in which either $ or A in 
(P^) is on shell: 



soft regime: 


As-. 


ko 


~ \k\ ~ 


Mv , 




potential regime: 




T 




\p\ ~ Mv , 


(2.4) 


ultrasoft regime: 


An-. 


ko 


~ \k\ ~ 


Mf2 





Ultrasoft gluons A^ are emitted as bremsstrahlung or from excited states in the bound 
system. Soft gluons As do not describe bremsstrahlung: Because in- and outgoing quarks 
$p are close to their mass shell, they have an energy of order Mv'^. Therefore, overall energy 
conservation forbids all processes with outgoing soft gluons but without ingoing ones, and 
vice versa, as their energy is of order Mv. 

The list of particles is not yet complete: In a bound system, one needs gluons which 
change the quark momenta but keep them close to their mass shell: 

Ap : fco ~ Mv^ , \k\r^ Mv (2.5) 

So far, only potential gluons and quarks, and ultrasoft gluons had been identified in the 
literature of power counting in NRQFT ^, |^. That the soft regime was overlooked 
cast doubts on the completeness of NRQFT after Beneke and Smirnov |^ demonstrated 
its relevance near threshold in explicit one- and two-loop calculations. In this article, the 
fields representing a non-relativistic quark or gluon came naturally by identifying all possible 
particle poles in the non-relativistic propagators, given the two scales at hand. 

When a soft gluon As couples to a potential quark $p, the outgoing quark is far off its 
mass shell and carries energy and momentum of order Mv. Therefore, consistency requires 
the existence of quarks in the soft regime as well, 

$s : T \p\ Mv . (2.6) 

As the potential quark has a much smaller energy than either of the soft particles, it can - 
by the uncertainty relation - not resolve the precise time at which the soft quark emits or 
absorbs the soft gluon. So, we expect a "temporal" multipole expansion to be associated 
with this vertex. In general, the coupling between particles of different regimes will not 
be point-like but will contain multipole expansions for the particle belonging to the weaker 
kinematic regime. For the coupling of potential quarks to ultrasoft gluons, this has been 
observed in Refs. 0. 

Propagators will also be different from regime to regime: For soft quarks, ^ is negligible 
against the kinetic energy T, so that the soft quark propagator may be expanded in powers 
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of and $s is expected to become static to lowest order. As the energy of potential 
gluons is much smaller than their momentum, the Ap-propagator is expected to become 
instantaneous. 

With these five fields $s, •I'p, Ag, Ap, A^^ representing quarks and gluons in the three dif- 
ferent non-relativistic regimes, soft, potential and ultrasoft, NRQFT becomes self-consistent. 
The application of these ideas to NRQCD with the inclusion of fermions and gauge particles 
is straightforward and will be summarised in the next publication ||^. An ultrasoft quark 
(which would have a static propagator) is not relevant for this paper. It is hence not con- 
sidered, as is a fourth ("exceptional") regime in which momenta are of the order Mv"^ and 
energies of the order Mv or any regime in which one of the scales is set by M. They do 
not derive from poles in propagators, and hence will be relevant only under "exceptional" 
circumstances. A future publication has to prove that the particle content outlined is not 
only consistent but complete. 

It is worth noticing that the particles of the soft regime can neither be mimicked by 
potential gluon exchange, nor by contact terms generated by integrating out the ultraviolet 
modes: Fields in the soft regime have momenta of the same order as the momenta of the 
potential regime, but much higher energies. Therefore, seen from the potential scale they 
describe instantaneous but non-local interactions, as pointed out in 0. Integrating out 
the scale Mv, one arrives at soft gluons and quarks as point-like multi-quark interactions 
in the ultrasoft regime. The physics of potential quarks and gluons will still have to be 
described by spatially local, but non-instantaneous interactions. The resulting theory - 



baptised potential NRQCD by Pineda and Soto |]T0[ - can be derived from NRQCD as 
presented here by integrating out the fields $8, and Ap. There is no overlap between 
interactions and particles in different regimes. 

In order to clarify this point, and before investigating the interactions of the various regimes 
further, the following example will demonstrate the power of dimensional regularisation in 
NRQFT. It also highlights some points which simplify the discussion of the following sections. 
The integral corresponding to a one- dimensional loop 

J(a, b) := J dk ^2 _ ^2 + i, k'-b^ + ie = ab{a + b) ^^'^^ 

is easily calculated using contour integration. Assuming := |^ -C 1, the dominating 
contributions come from the regions where \k\ is close to a ("smaller regime") or b ("larger 
regime"). Then, one can approximate the integral by 



I{a,b) 



+ 

\k\^a |fc|~f) 



In the first integral, k is small against b, so that a Taylor expansion in | ~ f in that regime 
is applicable and yields 

E / — ^ . (2.9) 



|fc|~a 



If k'^ becomes comparable to 6^, the expansion breaks down, so that the approximated 
integral cannot be solved by contour integration. In general, the (arbitrary) borders of the 
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integration regimes (the "cutoffs") will enter in the result, and lead to divergences as they 
are taken to infinity because of contributions from regions where \k\ ^ b ^ a. A cutoff 
regularisation may hence jeopardise power counting in v. 

Dimensional regularisation overcomes this problem in a natural and elegant way: If one 
treats ( p.9| ) as a (i-dimensional integral with d ^ 1 only at the end of the calculation, the 
exact result will emerge as a power series in f = |. First, one extends the integration regime 
from the neighbourhood of \a\ to the whole (i-dimensional space. Then, one calculates the 
integral order by order in the expansion, still treating ~ f ^ as formally small. Rewriting 



,2n 




(2.10) 



only the (m = ?T,)-term contributes thanks to the fact that dimensionally regularised integrals 
vanish when no intrinsic scale is present. 



/ 



fc" = . (2.11) 



(2vr) 
The result, 

is exactly the contribution one obtains in the contour integration from the pole at \k\ = a. 
Albeit the integral was expanded over the whole space, dimensional regularisation missed 
the poles at ±6 after expansion. 

The integration about \k\ ~ 6 is treated likewise by expansion and term-by-term dimen- 
sional regularisation. Adding this contribution, 

-iTT ^ a^" / -ivr 



to (|2.12| ), one obtains term by term the Taylor expansion of the exact result (|2.7|) in the 
small parameter v = ^. Each of the two regularised integrals sees only the pole in either of 
the regimes \k\ ~ a or \k\ ~ 6. Indeed, the overlap of the two regimes is zero in dimensional 
regularisation, even for arbitrary v. But then, the expansion in the two different regimes can 
be terminated only at the cost of low accuracy. 

One could therefore have started with the definition of two different integration variables, 
one formally living in the smaller regime with \Ka\ ^ a ^ vb, the other formally living in 
the larger regime, \Kh\ ~ 6: 

"dk ^ f dX + / d% (2.14) 



The momentum k is represented in each of the kinematic regimes by either Ka or Kh- The 
integrands must then be expanded in a formal way as if \Ka\ ^ a ^ vb and \Kf,\ ~ b. 
Otherwise, the poles are double counted. If one wants to calculate to a certain order in 
V, the expansion in the different variables ^ and is just terminated at the appropriate 
order. No double counting of the poles can occur. Coming back to the three different regimes 
of NRQFT (|2.4|), there will therefore be no double counting between any pair of domains. 
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Finally, note that the limit a — ^> is not smooth: For a = 0, dimensional regularisation 
of ( |2.9| ) is zero because of the absence of a scale ( ^.111 ). A pinch singularity encountered 
in contour integration at k = ±ie behaves hence like a pole of second order in dimensional 
regularisation and is discarded, see also App. 

By induction, the arguments presented here can be extended to prove that for any con- 
vergent one-dimensional integral containing several scales, Laurent expansion about each 
saddle point and dimensional regularisation gives the same result as contour integration. A 
formal proof of the validity of threshold expansion does presently not exist for the case of 
multi-dimensional and divergent integrals, but Beneke and Smirnov [^] could reproduce the 
correct structures of known non-trivial two-loop integrals using threshold expansion, which 
is highly suggestive that such a proof can be given. This claim is supported by the observa- 
tion that threshold expansion is very similar to the asymptotic expansion of dimensionally 
regularised integrals in the limit of loop momenta going to infinity, for which such a proof 
exists Q 



3 Rescaling Rules, Lagrangean and Feynman Rules 

In order to establish explicit velocity power counting in the NRQFT Lagrangean, one rescales 
the space-time coordinates such that typical momenta in either regime become dimensionless, 
as first proposed in M for the potential regime, and in p| for the ultrasoft regime: 



soft 
potential 
ultrasoft 



t={Mv)-^T, , x={Mv)-^X, , 
t={Mv^)-^T^ , x={Mv)~^Xs , 
t={Mv^)~^T^ , x={Mv'^)-^X^ 



(3.1) 



In order for the propagator terms in the NRQFT Lagrangean to be properly normalised, one 
has to set for the representatives of the gluons in the three regimes 



soft 
potential 
ultrasoft 



A,{x,t) = (Mv) As{Xs,T,) , 
Ap{x,t) = (Mv'i) Ap{X„T^) 
A^{x, t) = {Mv^) A(X„ T,) 



(3.2) 



and for the quark representatives 

soft: ^s{.x,t) 
potential: $p(x, t) 



(Mt;)2 0,(X3,T,) 



The rescaled free Lagrangean in the three regions reads then 



soft: 
potential: 
ultrasoft: 



v^2\ . 1 



d^Xs dTs 
d?X^ dr„ 



d3X„dT„-(5^A)(5M„) 



+ -(9^A)(9Ms 



am indebted to M. Beneke for conversation on this point 



(3.3) 



(3.4) 
(3.5) 
(3.6) 
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Here, as in the following, the positions of the fields have been left out whenever they coincide 
with the variables of the volume element. Derivatives are to be taken with respect to the 
rescaled variables of the volume element. The (un-rescaled) propagators are 

{T,P) i 

soft: $s : ^= = , (3.7) 

As : wm/mA 



F + ie 



potential: $p : » = ^ — , (3.8) 



i 



Ap-. 

k i 

ultrasoft: Ay^ : 'v/^/v^>^v^>^^ = — . (3.9) 

+ le 



As expected, the soft quark becomes static, resembling the quark propagator in heavy quark 
effective theory, and the potential gluon becomes instantaneous. In order to maintain velocity 
power counting, corrections of order v or higher must be treated as insertions as in the 



example, ( |2.7|) . Insertions are represented by the (un-rescaled) Feynman rules 



k 



=^^= = -''f^ = ^(^) ' — -X— = +■< = 0{v') . (3.10) 

Except for the physical gluons A^ and Au, there is no distinction between Feynman and 
retarded propagators in NRQFT: Antiparticle propagation has been eliminated by the field 
transformation from the relativistic to the non-relativistic Lagrangean, and both propagators 
have maximal support for on-shell particles, the Feynman propagator outside the light cone 
vanishing like e~*^. Feynman's perturbation theory becomes more convenient than the time- 
ordered formalism, as less diagrams have to be calculated . 

Finally, the interaction part of the Lagrangean reads (neglecting for the moment the $^ 
vertex in (p^) ) 

soft: d'Xs dTs i-g) [(A + A(^s, vT^) + v A^ivX^, vTs))4(f)s + (3.11) 

+ (A0sVp(^s, vTs) + h.c] 

3-r AT (-n\f—A J- A („Y T\\A.t 



potential: d^X^ dT, {-g)(—Ap + A^ivX^, r,) UUp . (3.12) 



Note that the scaling regime of the volume element is set by the particle with the highest 
momentum and energy. Vertices like ^s0p0p cannot occur as energy and momentum must 
be conserved within each regime to the order in v one works. Amongst the fields introduced, 
these are the only interactions within and between different regimes allowed. One sees that 
technically, the multipole expansion comes from the different scaling of x and t in the three 
regimes. It is also interesting to note that there is no choice but to assign one and the 
same coupling strength g to each interaction. Different couplings for one vertex in different 
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regimes are not allowed. This is to be expected, as the example ( |2.14|) demonstrated that the 
fields in the various regimes are representatives of one and the same non-relativistic particle, 
whose interactions are fixed by the non-relativistic Lagrangean ( p.2|) . 
The interaction Feynman rules are 



(T,P) {T',p') 



]k = -igi2n)UiT + r + ko)5^^\p+p' + k) = 0{-. 



(3.13) 



d 



i(7(27r)4 5(T + T' + A;o)[exp(fc-^^^-^) + = o{^^) ,(3.14) 




-i^(27r)^(5(T + T' + A;o) +^) = 0(w°) , 



(3.15) 



-i^(27r)^ 



exp 



(t' 



d 



d{T + ko 



-ig{2TTy exp (jt, 



d 



d{T + T') 



5{T + k^)\ S^^\p+p' + k) = C(e'') ,(3.16) 



5(T + T')1 6^^\p+p' + k) = 0(v^eX^.17) 



exp fco 



{k 



d 



X 



d{T + T 
exp(£ ' 



5{T + T') 



X 



(3.18) 



d{p + p'] 



6^^\p+p')\ = Give") 



The exponents representing the multipole expansion have to be expanded to the desired order 
in V. Double counting is prevented by the fact that in addition to most of the propagators, 
all vertices are distinct because of different multipole expansions. 

Using the equations of motion, the temporal multipole expansion may be re- written such 
that energy becomes conserved at the vertex. Now, both soft and potential or ultrasoft 
energies are present in the propagators, making it necessary to expand it in ultrasoft and 
potential energies. An example would be to restate the vertex ( p.lTj) as 



-ig{2TT)UiT + r + kp^o)5^^\p+p' + k) = OiV^) 
I 

and the soft propagator as containing insertions 0{v) for potential energies kp 



(3.19) 



(T+fcp,o,p) 



T + ie 



E 

n=0 



T 



(3.20) 



The same holds of course for the momentum-non-conserving vertices. 

In the renormalisation group approach, there is therefore only one relevant coupling (i.e. 
only one which dominates at zero velocity): As expected, it is the $p$py4p coupling providing 
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the binding. The "I's^s^p-; ^'s'^'s^u-couphngs and both insertions ( p.lO| ) are irrelevant. The 
marginal couplings $p$pAu, $s'^'s^s and $s'^'p^s are irrelevant in gauge theories in carefully 
chosen gauges like the Coulomb gauge. This point will be elaborated upon in the future 0. 

The velocity power counting is not yet complete. As one sees from the volume element 
used in ( |3.11|) , the vertex rules for the soft regime count powers of v with respect to the 
soft regime. One hence retrieves the velocity power counting of Heavy Quark Effective 
Theory |T2| , |13[ (HQET), in which the interactions between one heavy (and hence static) 
and one or several light quarks are described. Usually, HQET counts inverse powers of mass 
in the Lagrangean, but because in the soft regime Mv ~ const., the two approaches are 
actually equivalent. HQET becomes a sub-set of NRQCD, complemented by interactions 
between soft (HQET) and potential or ultrasoft particles. 

In NRQCD with two potential quarks as initial and final states, the soft regime can occur 
only inside loops, as noted above. Therefore, the power counting in the soft sub-graph has to 
be transfered to the potential regime. Because soft loop momenta scale like [d^/cg] ~ v^, while 
potential ones like [d^A;p] ~ f^, each largest sub-graph which contains only soft quarks and 
no potential ones (a "soft blob") is enhanced by an additional factor -. Couplings between 
soft quarks and any gluons inside a blob take place in the soft regime and hence are counted 
according to the rules of that regime. Each soft blob contributes at least four orders of 
but only one inverse power of f ~ g^. Power counting is preserved. These velocity power 
counting rules in loops are verified in explicit calculations of the exemplary graphs (see also 
below), but a rigorous derivation is left for a future publication [§]. 

With rescaling, multipole expansion and loop counting, the velocity power counting rules 
are established, and one can now proceed to check the validity of the proposed Lagrangean, 
matching NRQFT to the relativistic theory in the examples given by Beneke and Smirnov 



4 Model Calculations 

The first example is the lowest order correction to the two quark production graph. With- 
out proof, it will be used that in dimensional regular isat ion, one can match NRQFT and 
the relativistic theory graph by graph, so that not the whole scattering amplitude has to 
be considered [^. The collection of graphs to be matched to the relativistic diagram is 
depicted in Fig. |1|. Here and in the following, hard (ultraviolet) contributions will not be 




Figure 1: Matching of the 0{g'^) correction to two quark production oS an external current 
to lowest order in v in each of the three regimes. 

shown explicitly. They are taken care of by the four-quark interaction of the non-relativistic 
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Lagrangean ( |2.2| ) and renormalisation of the external currents . 

The energy and momentum routing has been chosen to be the one of the non-relativistic 
center of mass system, with 2T the total kinetic energy, and y = — (p)^ oc —v^ the relative 
four-momentum squared of the outgoing quarks as indicator for the thresholdness of the 
process considered. Thanks again to dimensional regularisation, any other assignment can 
be chosen and reproduces the result. 

The vanishing of the ultrasoft gluon exchange diagram and the value of the potential 
gluon exchange diagram have already been calculated in 0. The soft exchange diagram 
vanishes, so that no new contribution is obtained. It is not even necessary to specify how 
soft quarks couple to external sources: If energy is conserved at the production vertex, the 
integral to be calculated is 



1 



1 



(4.1) 



As the gluon is soft, T <^ k^ and the quark propagators must be expanded in T/k^ ~ u, 
giving zero to any order as no scale is present in the dimensionally regularised integral. If 
energy is not conserved at the production vertex, the soft quark propagator is and the 
contribution vanishes again. Therefore, there is no coupling of soft subgraphs to external 
sources to any order in v. Soft quarks in external lines are far off their mass shell and 
hence violate the assumptions underlying threshold expansion and NRQFT. In general, we 
conclude that soft quarks are present only in internal lines, and that the first non- vanishing 
contribution from the soft regime for the production vertex occurs not earlier than at 0{g^). 

The first soft non-zero contribution comes actually from the two gluon direct exchange 
diagram of Fig. |] calculated by Beneke and Smirnov using threshold expansion. The 
Mandelstam variable t = —{p — p'Y describes the momentum transfer in the center of mass 
system. The ultraviolet behaviour of this graph is mimicked in NRQFT by a four-fermion 



(T,p) 




{T,p') 


c 







fcT c 




\ 1 {ka, ^ 
1 p—'p' +k) 


o 
o 









{T,-P) 



iT,-p') 



(fco,p+fc) 




(-feo,-p-fc) 

Figure 2: Planar 0{g'^) contributions to Coulomb scattering. The four-point interaction and 
insertion diagrams are not displayed. 



exchange given by the vertex ic2 = 24-J^m'2 ~ of the Lagrangean 

the rescaling rules is seen to be 0{v). 



which using 
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The Feynman rules (|3.14| ) give that the y4uv4u-diagram is of order e^ with a leading loop 
integral contribution (similar to fi. (32)]) 

d* 1 1 1 1 



The diagram is expected to be zero to all orders since the ultrasoft gluons do not change 
the quark momenta and therefore the scattering takes place only in the forward direction, 
p = p. Upon employing the on-shell condition for potential quarks, T = ^ to leading order, 

it indeed vanishes as no scale is present. Since T — ^ ~ Mv'^ <^ |A;| ~ Mv (and ko ~ Mv'^) 
in the potential regime, this is a legitimate expansion. The AuAp and ApAu contributions 
{0{-e^)) are zero for the same reason. The lowest order contribution to the ApAp graph 

(Oih) is 

^ ^ 1 1 (4.3) 



{27rY p-ie (p-p' + ^)2_ie t + ko - + ie T-k^-^^ + ie 



In the light of the discussion at the end of Sect. H, it is most consistent to perform the k^ 
integration by dimensional regularisation, using / = J (L)d-a , cr ^ 1 0, Chap. 

4.1]. Split dimensional regularisation was introduced by Leibbrandt and Williams to cure 
the problems arising from pinch singularities in non-covariant gauges. Appendix ^ shows 
that in the case at hand, it has the same effect as closing the fco-contour and picking the quark 
propagator poles prior to using dimensional regularisation in c? — 1 Euclidean dimensions. To 
achieve 0{v^) accuracy, one must also consider one insertion ( p.lO|) at the potential gluon 
lines, giving rise to a contribution 



J— \ kl ( J— + \ ) X 

(27iY - ie (p _ + ky - ie ^ ^ P - ie {p - p' + k^ - ie ' 

X K — . (4.4) 

T + ko-^ + ieT-ko-^ + ie 



The fco integration is naively linearly divergent, and hence closing the contour is not straight- 
forward. As App. ^ demonstrates, split dimensional regularisation circumvents this problem. 
The sum of both contributions ([4.3|/[4.4|) , 



7E-ln- — - , (4.5) 



agrees with fl. (31)] when one keeps in mind that in that reference, heavy particle external 
lines were normalised relativistically, while a non-relativistic normalisation was chosen here. 
Also, this article uses the MS rather than the MS scheme. Near threshold, the scale is set 
by the total threshold energy Air/i'^ = 4(M + T)^. 

The soft gluon part is to lowest order (C(f ~^) because of one soft blob) given by 

d'^k 1 1 11 , 

(4.6) 



{27rY kl-k^ + ie k^-{p-p' + kf + ie + ie -^o + ie 
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which corresponds to p, fl. (33)]. Now, spht dimensional regularisation must be used if 
no ad-hoc prescription for the pinch singularity at /cq = is to be invoked. That the 
pinch is accounted for by potential gluon exchange and hence must be discarded, agrees 
with the intuitive argument that zero four-momentum scattering in QED is mediated by a 
potential only, and no retardation or radiation effects occur. On the other hand, the model 
Lagrangean contains three marginal couplings as seen at the end of Sect. |^, which may give 
finite contributions as energies and momenta of the scattered particles go to zero. Although 
the prescription and the result from split dimensional regularisation coincide in the present 
case as demonstrated at the end of Sect. ^ and in App. this may not hold in general. The 
result to 0{v^) exhibits another coUinear divergence. 



7e 



In 



+ 



247r2M2 



1 + 



2y 



t M-rf 



- 7e - In 



47r/i^ 



(4.7) 



and agrees with the one given by Beneke and Smirnov 0, fi. (36)]. The second term comes 
from insertions and multipole expansions to achieve 0{v) accuracy. 

It is easy to see that the power counting proposed works. As expected, the potential 
diagram is ^/y oc v stronger that the leading soft contribution, and t^/y oc stronger than 
the four-fermion interaction. 

In conclusion, the proposed NRQFT Lagrangean reproduces the result for the planar 
graph of the relativistic theory only if the soft gluon and the soft quark are accounted for: 
The four-fermion contact interaction produces just a -ga-term, graphs containing ultrasoft 
gluons were absent, and the potential gluon ( [4.5[ ) gave no 0{y^) contribution. This shows 
the necessity of soft quarks and gluons. The coupling strength of the $s^s'^'p vertex is also 
seen to be identical to the other vertex coupling strengths, g. 

The planar fourth order correction to two quark production (Fig. ^ was also compared 
to the result of [^, and is correctly accounted for when the Feynman rules proposed above 
are used to 0{v^). 

{T,P) 



(2T,0) 





+ 




-P) 



Figure 3: The non-vanishing contributions to the planar fourth order correction to two quark 
production. Diagrams with insertions or four-point interactions not displayed. 



5 Conclusions and Outlook 

The objective of this article was a simple presentation of the ideas behind explicit velocity 
power counting in dimensionally regularised NRQFT. It started with the identification of 
three different regimes of scale for on-shell particles in NRQFT from the poles in the non- 
relativistic propagators. This leads in a natural way to the existence of a new quark field 
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and a new gluon field in the soft scaling regime ~ \p\ ~ Mv. In it, quarks are static 
and gluons on shell, and HQET becomes a sub-set of NRQCD. Neither of the five fields in 
the three regimes should be thought of as "physical particles". Rather, they represent the 
"true" quark and gluon in the respective regimes as the infrared-relevant degrees of freedom. 
None of the regimes overlap. An NRQFT Lagrangean has been proposed which leads to the 
correct behaviour of scattering and production amplitudes. It establishes explicit velocity 
power counting which is preserved to all orders in perturbation theory. The reason for the 
existence of such a Lagrangean, once dimensional regularisation is chosen to complete the 
theory, was elaborated upon in a simple example: the non-commutativity of the expansion 
in small parameters with dimensionally regularised integrals. 

Due to the similarity between the calculation of the examples in the work presented here 
and in 0, one may get the impression that the Lagrangean presented is only a simple re- 
formulation of the threshold expansion. Partially, this is true, and a future publication 
will indeed show the equivalence of the two approaches to all orders in the threshold and 
coupling expansion. A list of other topics to be addressed there contains: the straightforward 
generalisation to NRQCD; a proof whether the particle content outlined above is not only 
consistent but complete, i.e. that no new fields (e.g. an ultrasoft quark) or "exceptional" 
regimes arise; an investigation of the influence of soft quarks and gluons on bound state 
calculations in NRQED and NRQCD; a full list of the various couplings between the different 
regimes and an exploitation of their relevance for physical processes. The formal reason why 
double counting between different regimes and especially between soft and ultrasoft gluons 
does not occur, a derivation of the way soft quarks couple to external sources, and the role 
of soft gluons in Compton scattering deserve further attention, too. 

I would like to stress that the diagrammatic threshold expansion derived here allows 
for a more automatic and intuitive approach and makes it easier to determine the order in 
y/—y (X V to which a certain graph contributes. On the other hand, the NRQFT Lagrangean 
can easily be applied to bound state problems. As the threshold expansion of Beneke and 
Smirnov starts in a relativistic setting, it may formally be harder to treat bound states there. 
Indeed, I believe that even if one may not be able to prove the conjectures of the one starting 
from the other, both approaches will profit from each other in the wedlock of NRQFT and 
threshold expansion. 
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A Some Details on Split Dimensional Regularisation 

This appendix presents the part of the calculations in the examples of Sect. |^ which makes 
use of split dimensional regularisation as introduced by Leibbrandt and Williams 0. In its 
results, split dimensional regularisation agrees with other methods to compute loop integrals 



in non-covariant gauges, such as the non-principal value prescription |15|, but two features 
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make it especially attractive: It treats the temporal and spatial components of the loop 
integrations on an equal footing, and no recipes are necessary. Rather, it uses the fact that, 
like in ordinary integration, the axioms of dimensional regularisation [|1^, Chap. 4.1] allow 
to split the integration into two separate integrals: 

d'% /• d% d'^-^k 

(A.l) 



{2ttY J {2TTy {2tt 



\d—(T 



Both integrations can be performed consecutively, and the limit cr — 1 can - if finite - be 
taken immediately, because the integration over the spatial components of the loop momen- 
tum in ([A.l|) is still regularised in d — 1 dimensions. Finally, the limit ci — > 4 is taken at the 
end of the calculation. 

Equation (^73|) contains the simplest ko sub-integral: 



d% 1 1 



(A.2) 



Using standard formulae for dimensional regularisation in Euclidean space [|16| , App. B], the 
result is finite as o" — 1: 



/ 



d% 1 _ r[l-|] / jk + p)^ . 

(27r)- [T-^ + ie]2 - (4vr)-/2r[l] I L 2M 

4(^-^^ + i^)" asa-1 (A.3) 



2 V 2M 



It is no surprise that closing the contour produces the same result, because for any finite 
integral, the answer of all regularisation methods have to coincide. The integral over the 
spatial components of the loop momentum is now straightforward. 

The potential gluon diagram with one insertion at a gluon leg (|4.4| ) yields a split di- 
mensional integral which diverges linearly in ko, so that naive contour integration is not 
legitimate. 

^0 . ' fT-^l±^) as.^1 (A.4) 



(2^). lT-^ + ie]^-ki 2 



2M 



To arrive at this result, the numerator was re- written as (feg — (T — ^^j^Y) + (T — ■^^^^)^- 
Its first term cancels the denominator, yielding an integral without scale which therefore 
vanishes in dimensional regularisation. The second term has been calculated in ( |A.3D . The 
integral over the spatial components of the loop momentum provides again no complications, 
leading to ([4.5[). 

Finally, it was already shown at the end of Sect. |^ that dimensional regularisation discards 
pinch singularities encountered in contour integrations. This is validated again by looking 
at the split dimensional integral for k^ in the soft gluon contribution ( |4.6| ), 

' ' ' (A.5) 



{2ttY kl -a? kl- 62 -kl 
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where c? := fc^ — ie and 6^ := {p — p' + — ie. After combining denominators, the resulting 
integral is simple: 

CTj ^ -—3 i 1 11 

This agrees with the result of Beneke and Smirnov fi. (34)] who use contour integration and 
drop the contribution from the pinch singularity. The integral over the spatial components 
of the loop momentum provides again no unfamiliar complications, leading to ([4.7|) . 
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